Effects of electron-phonon coupling in the Kondo regime of a two-orbital molecule 
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We study the interplay between strong electron-electron and electron-phonon interactions within a two-orbital 
molecule coupled to metallic leads, taking into account Holstein-like coupling of phonons to the molecular 
occupancies as well as phonon-mediated interorbital tunneling. By combining canonical transformations with 
numerical renormalization-group calculations to address the interactions nonperturbatively and on equal footing, 
we obtain a comprehensive description of the system's many-body physics. The electron-phonon interactions 
strongly modify the bare orbital energies and the Coulomb repulsion between electrons in the molecule, and 
tend to inhibit tunneling of electrons between the molecule and the leads. The consequences of these effects are 
considerably more pronounced when both molecular orbitals lie near the Fermi energy of the leads than when 
only one orbital is active. In situations where a local moment forms on the molecule, there is a crossover with 
increasing electron-phonon coupling from a regime of collective Kondo screening of the moment to a limit of 
local phonon quenching. At low temperatures, this crossover is associated with a rapid increase in the electronic 
occupancy of the molecule as well as a marked drop in the linear electrical conductance. 

PACS numbers: 71.38.-k, 72.15.Qm, 72.10.Fk, 73.23.-b, 73.23.Hk, 73.63.-b, 73.63.Kv, 73.63.Rt 



I. INTRODUCTION 

Single-molecule junctions^ are structures consisting of a 
single molecule bridging the gap between source and drain 
electrodes, allowing electronic transport when a bias voltage 
is applied across the structure. These systems, which manifest 
a rich variety of experimentally accessible physics in a rela- 
tively simple setting,- have attracted much theoretical and ex- 
perimental effort in connection with molecular electronics. 6 ' 7 
A major goal of these efforts has been to take advantage of 
natural or artificial molecules for technological purposes. Ex- 
amples of single-molecule junctions encompass, for exam- 
ple, single hydrogen molecules 8 - 1 - and more complex struc- 
tures such as 4, 4'-bipyridine molecules coupled to metallic 
nanocontacts J^dl 

An important ingredient in transport through molecular 
systems is the electron-electron interaction (Coulomb repul- 
sion), the effect of which is greatly enhanced by the spa- 
tial confinement of electrons in molecules^ Electron-electron 
(e-e) interactions are known to produce Coulomb blockade 
phenomena 15 and Kondo correlations at low temperatures ji&ii 
Confined electrons are also known to couple to quantized vi- 
brations (phonons) of the molecules J 8 resulting in important 
effects on electronic transport.-^— Single-molecule junctions 
therefore provide a valuable opportunity for studying charge 
transfer in systems with strong competing interactions. 

It has recently been demonstrated that the energies of the 
molecular orbitals in a single-molecule junction can be tuned 
relative to the Fermi energy of the electrodes by varying 
the voltage applied to a capacitively coupled gate.— Simi- 
lar control has for some time been available in another class 
of nanoelectronic device: a quantum dot coupled to a two- 
dimensional electron gas. 29,30 The electrons confined in a 
quantum dot couple — in most cases quite weakly — to col- 
lective vibrations of the dot and its substrate.— In single- 
molecule devices, by contrast, the confined electrons interact 



with local vibration modes of the molecule that can produce 
pronounced changes in the molecular electronic orbitals. Con- 
sequently, electron-phonon (e-ph) interactions are expected to 
play a much more important role in molecules than in quan- 
tum dots. 

From the theoretical point of view, addressing both e-e and 
e-ph interactions from first principles is a very complicated 
task. However, simple effective models can provide good 
qualitative results, depending on the parameter regime and 
the method employed to solve the model Hamiltonian.— For 



example, the essential physics of certain experiment 



,21.23 



ap- 



pears to be described by variants of the Anderson-Holstein 
model, which augments the Anderson modelM for a magnetic 
impurity in a metallic host with a Holstein coupling 3 -^ of the 
impurity charge to a local phonon mode. The model has been 
studied since the 1970s in connection with other problems^^. 
prior to its application to single-molecule devices.^ 5 -^ 6 . Vari- 
ous analytical approximations as well as nonperturbative nu- 
merical renormalization-group calculations have shown that 
in equilibrium, the Holstein coupling reduces the Coulomb 
repulsion between two electrons in the impurity level, even 
yielding effective e-e attraction for sufficiently strong e-ph 
coupling. Increasing the e-ph coupling from zero can produce 
a smooth crossover from a conventional Kondo effect, involv- 
ing conduction-band screening of the impurity spin degree of 
freedom, to a charge Kondo effect in which it is the impurity 
"isospin" or deviation from half-filling that is quenched by the 
conduction band. 

Single-molecule devices at finite bias are usually studied 
via nonequilibrium Keldysh Green's functions that systemati- 
cally incorporate the many-body interactions within a system. 
Although this approach has proved to be the most reliable for 
calculation of transport properties, the results are highly sen- 
sitive to the approximations made. For instance, the equation- 
of-motion technique^ 7 - generates a hierarchy of coupled equa- 
tions for Green's functions containing 2« fermionic operators 



for n = 1,2,3,...: a hierarchy that must be decoupled at some 
level in order to become useful. The most commonly used de- 
coupling scheme is based on a mean-field decomposition of 
the n — 2 Green's functions, leading to the well-known Hub- 
bard I approximation. 4? This approximation give reasonable 
results for temperatures T above the system's Kondo temper- 
ature Tk, but as it neglects spin correlations between localized 
and conduction electrons, it fails in the Kondo regime. 

A few years ago, two of us applied the equation-of-motion 
method decoupled at level n — 2 to study a single-molecule 
junction that features phonon-assisted interorbital tunneling.— 
However, to capture the physics at T < Tk requires extension 
of the equation-of-motion hierarchy to higher order, which in 
most cases is carried out in the limit of infinite Coulomb in- 
teraction. The Kondo regime may also be studied via dia- 
grammatic expansion within the non-crossing approximation, 
which is again most straight forward in the infinite-interaction 
limits 

This paper reports the results of an investigation of the 
Kondo regime of a two-orbital molecule, with focus on sit- 
uations in which Coulomb interactions are strong but finite. 
In order to treat e-e and e-ph interactions on an equal ba- 
sis, we employ Wilson's numerical renormalization-group 
approach^— which allows complete access to the equilib- 
rium behavior and linear response of the system for tem- 
peratures all the way to absolute zero. We show that the 
renormalization of e-e interactions is strongly dependent on 
the energy difference between the two molecular levels. For 
small interorbital energy differences, the renormalization is 
significantly enhanced compared with the situation of one ac- 
tive molecular orbital considered in previous work. This en- 
hancement is detrimental for realization of the Kondo effect 
but improves the prospects for accessing a phonon-dominated 
regime of effective e-e attraction. A sharp crossover be- 
tween Kondo to phonon-dominated regimes has signatures in 
thermodynamic properties and in charge transport across the 
molecule. 

We note that we leave for future study cases involving two 
degenerate (or nearly degenerate) molecular orbitals lying be- 
low the chemical potential of the leads. In such cases, e-e 
interactions favor the presence of an unpaired electron in each 
orbital, and electron-phonon interactions may be expected to 
significantly affect the competition between total-spin-singlet 
and triplet ground states^— 

The rest of the paper is organized as follows: Section [XT] 
describes our model system and provides a preliminary analy- 
sis via canonical transformations. Section [HI] summarizes the 
numerical solution method and Sec. [IV] presents the results. 
Concluding remarks appear in Sec. \V\ 



H. MODEL AND PRELIMINARY ANALYSIS 

A. Model Hamiltonian 

We consider a system composed of a two-level molecule 
interacting with a local phonon mode and also coupled to two 
metallic leads, as shown schematically in Fig.Q] This system 
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FIG. 1 : (Color online) Schematic representation of the model studied 
in this work. A molecule with two active orbitals (a and P) spans the 
gap between left and right electrodes. The molecular orbitals are 
subject both to e-e and e-ph interactions. 

is modeled by the Anderson-type Hamiltonian 

H = H mo ] + //leads + Z/mol-leads> (1) 

with H mo \ describing the isolated molecule, Z/i ea d s modeling 
the leads, and Z/ mo i-ieads accounting for electron tunneling be- 
tween the molecule and the leads. 

The molecular Hamiltonian can in turn be divided into four 
parts: H mo \ — H e + Z/ p h + //hoi + //tun- Here, the electronic part 

He = ^ (e,«; + f/,«, T «,u) + U'n a np, (2) 

i=afi 

where n !(r = dj d. is the number operator for electrons of 
energy e, and spin <x in molecular orbital i — a or f3, rii = + 
riii, U, and U' parametrize intraorbital and interorbital 
Coulomb repulsion, respectively. Without loss of generality, 
we take sp > s a . The phonon part 

Z/ph = two n b (3) 

describes a dispersionless optical phonon mode of energy 
Hcjq, with «i = tfb. The remaining two parts of H mo \ describe 
two different types of e-ph interaction: 

//Hoi = MnolC^ + b) (4) 

is a Holstein coupling between the phonon displacement and 
the combined occupancy (i.e., charge) 

«mol = n a + lip (5) 

of the two molecular orbitals, while 

Htun = A' YjdLdpa + ^Ar)^ + b) (6) 
cr 

describes interorbital tunneling accompanied by emission or 
absorption of a phonon. 

The left {( = L) and right (I = R) leads are represented by 

//leads = J] £k4 ka C (kiT , (7) 

£,k,cr 

where C( k(r annihilates an electron with energy s k , wave vector 
k, and spin cr in lead I. For simplicity, each lead is character- 
ized by a flat density of states 
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where N s is the number of lattice sites in each lead, D is the 
half bandwidth and 0(x) is the Heaviside function. 
Lastly, 

ff m0 l-leads = Z Z M^to + V«4Ar) (9) 

i=arjS S,k,<r 

describes tunneling between the molecular orbitals and the 
leads, allowing transport through the system. We assume that 
the tunneling matrix elements are real and have left-right sym- 
metry so we can write Vu — Vtf = V,. Then it is useful to 
perform an even/odd transformation 



C e ko- - —p(cRkcr + CLktr) 
V2 

Coko- = —^(cRka- ~ Cikcr), 

V2 

which allows Eq. (O to be rewritten 



(10) 
(ID 



mol- leads 



eko~ 



■ c\ d. 

ewer tcr 



), (12) 



k.<r 



With this transformation, the odd-parity degrees of freedom 
are fully decoupled from the molecular orbitals, and can safely 
be dropped. As a result, the problem reduces to one effective 
conduction channel described by a modified 



=Z- 



(13) 



This channel is still described by the density of states in Eq. 
(0, and it imparts to molecular orbital i a width 



F, = nVf/D. 



(14) 



A similar transformation to an effective one-channel model 
can be derived in any situation where the tunneling matrix el- 
ements satisfy VlcVr/s = Vl^Vro, ensuring that both molecu- 
lar orbitals couple to the same linear combination of left- and 
right-lead states. 

Since this model of a two-orbital single-molecule junction 
contains eleven bare energy scales, we will focus our nu- 
merical calculations on a few examples that illustrate inter- 
esting physics, rather than attempting a complete exploration 
of the model's parameter space. It should be noted that in 
the limit where one of the molecular orbitals (B, say) is re- 
moved or becomes decoupled from the rest of the system, the 
model reduces to the Anderson-Holstein Hamiltonian studied 
previously^i^ 



B. Preliminary analysis via canonical transformation 

Insight can be gained into the properties of the two-orbital 
model by performing a canonical transformation of the Lang- 
Firsov type££ from the original Hamiltonian ([T) to H = 



e Sl He Sl , Following extensive algebra, it can be shown that 
the choice 



S 1 = T — «moi(b 1 - b) 



(15) 



eliminates the Holstein coupling between the local phonons 
and the molecular electron occupancy [Eq. ©], leaving a 
transformed Hamiltonian 



H = H e + H p h + H mn + 7/leads + ffmol-leads, 



(16) 



in which 7/ p h and //i ea ds remain as given in Eqs. (01 and ( fT3l , 
respectively; H e has the same form as H e [Eq. (f2|] with the 
replacements 



A /tiajQ, 



Ui -> Ui = Ui - 2A 2 lhu Q , 
U' -> U> = V - 2A 2 /hhj Q ; 

the interorbital tunneling maps to 



(17a) 
(17b) 
(17c) 



where cr 
comes 



b r +b- ^-(1 + n ad . + tips) 
(18) 

-<t\ and the molecule-leads coupling term be- 



mol -leads 



— V Vi y {B\ct. c . +B 3 c\ d. ), (19) 



with 



B s = exp 



ttiOQ 



(b ] ' - b) 



= Bl 



(20) 



This transformation extends the one applied previously 
(e.g., see, Ref. l39l) to the Anderson-Holstein model. It ef- 
fectively eliminates the Holstein Hamiltonian term [Eq. (|4]i] 
by mapping the local phonon mode to a different displaced 
oscillator basis for each value of the total molecular occu- 
pancy « mo i, namely, the basis that minimizes the ground-state 
energy of H e + H p h + Hu i- There are three compensating 
changes to the Hamiltonian: (1) Shifts in the orbital energies 
[Eq. ( |17at l and, more notably, a reduction in the magnitude — 
or even a change in the sign — of each e-e interaction within 
the molecule [Eqs. (117bb and d!7cH . These renormalizations 
reflect the fact that the Holstein coupling lowers the energy 
of doubly occupied molecular orbitals relative to singly oc- 
cupied and empty orbitals. This well-known effect under- 
lies the standard e-ph mechanism for superconductivity. (2) 
Addition of correlated (molecular-occupation-dependent) in- 
terorbital tunneling [the /l-dependent term in Eq. ( [181 )1 to the 
phonon-assisted tunneling present in the original Hamilto- 
nian. (3) Incorporation into the molecule-leads coupling [Eq. 
( fT9l )l of operators B , and B^ that cause each electron tunnel- 
ing event to be accompanied by the creation and absorption of 
a packet of phonons as the local bosonic mode adjusts to the 
change in the total molecular occupancy n mo \. 
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The effects of the phonon-assisted interorbital tunneling 
term // tun can be qualitatively understood by rewriting Eq. 
(fT6b in terms of even and odd linear combinations of the a 
and B molecular orbitals: 



deer - —=(d aC r + dpo), d ocr - — — (d acr - dpa). (21) 

V2 V2 



In this parity basis, Eq. ( fTST l becomes 



#tun = A' ^J,n m - n oa ) 



t , 2A . 
V +b - - — (1 + n ed . + n oB ) 
TiaiQ 



, (22) 



where n pa = d pa -d piT for p — e or o. The phonon-assisted 
tunneling component of H tun (i.e., the original // tun ) can be 
eliminated by performing a second Lang-Firsov transforma- 
tion 



H = e^He 



S 1 U„- S 2 



(23) 



with 



S 2 = - — (n e - n„)(b' 1 ' - b), 
ncoo 



(24) 



where n p - n p y + n p ±. Lengthy algebra reveals a transformed 
Hamiltonian 



H - H e + H p h + //leads + //mol-leads, 



(25) 



where 



H e = X! i^P n P + Up n pt n pl) + J^j&W 

fleer W oo~ U±_ fteo~fto<rj 

p=e,o or 

+ B 2A' d lad e a) ( 26 ) 

cr 

+ j(s;s- +s;s; - B\ v t e r a - B Av r e t ) 



with S* = (S;) f = c^c pl and I* = (I p )"'" = c^c^ being spin- 
and charge-raising operators, respectively, and 



//mol-leads - ~= ^ Vp X^L^jLr^kcr + B ^' C LrV)- 
V" J p=e,o k,tr 

(27) 

The renormalized parameters and operators entering Eqs. 



and ( |27] | are 



2 Swo ' 
2C7' + £/ a + Up 2A 2 p 
TiujQ ' 



#11 = ' ~ 



4 

2A P A n 



2U' + U a + Up 2A e A a 
4 ttu>o 



I 

W = 
J = 
V = 



2 




U a - 


Up 


2 




2U' - 


-Ua 
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V a ± 


Vp 



where 



A eo — A ± A . 



(28a) 

(28b) 
(28c) 
(28d) 
(28e) 
(28f) 
(28g) 
(28h) 

(29) 



Since the e-e interactions are expressed much less compactly 
in the parity basis than in the original basis of a and B orbitals, 
the elimination of the boson-assisted interorbital tunneling 
from the Hamiltonian comes at the price of much greater com- 
plexity in H e compared to H e [Eq. d2J] and H e . It is no- 
table, though, that the e-e repulsion between two electrons 
within the even-parity [odd-parity] molecular orbital under- 
goes a non-negative reduction proportional to A% = (A + A') 2 
[A 2 , — (A - A') 2 ]. By contrast, the Coulomb repulsion be- 
tween electrons in orbitals of different parity undergoes a shift 
proportional to —A e A a = A' 2 - A 2 that may be of either sign. 
Whereas large values of |A| favor double occupancy of both 
the a and the B molecular orbital, large values of \A'\ favor 
double occupancy of either the e or the o linear combination 
[the degeneracy between these alternatives being broken by an 
amount (2s e + U e ) - (2s + U ) = -16AA' /%a>o\- Both lim- 
its yield a unique many-body ground state of a very different 
character than the spin-singlet Kondo state. 

Since S\ defined in Eq. dT3T > can be rewritten Si = 
(A/ha>o)(n e + « )(£>' - b), it commutes with S2 given in Eq. 
( l24b . As a result, the two Lang-Firsov transformations can be 
combined into a single canonical transformation 



H = e s He~ s 



with 



S = Si +s 2 



halo 



'-(b 1 ' - b). 



(30) 



(31) 



This canonical transformation maps the original phonon anni- 
hilation operator b to 



f St -S 1 ^lle + A H 

b = e be =b- 



(32) 
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If the phonon energy Tiojq were to greatly exceed the thermal 
energy k B T and all other energy scales within the model, the 
system's ground state would be expected to be characterized 
by b ' b = and, hence, 



(n b ) = (tfb) = 



lt A e n e + /t » \ 2 \ 



(33) 



However, the physical limit of greater interest is one in which 
the Coulomb interactions U a , Up, and U' are larger than fru)Q. 
The applicability of Eq. (l33l to such situations is addressed in 
Sec. |IV] 



III. NUMERICAL RENORMALIZATION-GROUP 
APPROACH 

In order to obtain a robust description of the many-body 
physics of the model, we treat the Hamiltonian (Q~|i using Wil- 
son's numerical renormalization-group (NRG) method*^— 
as extended to incorporate local bosonic degrees of freedom. 41 
The effective conduction band formed by the even-parity com- 
bination of left- and right-lead electrons is divided into loga- 
rithmic bins spanning the energy ranges DA~ ( "' +I) < ±e < 
DA '" for m - 1, 2, 3, . . ., for some discretization parameter 
A > 1 . After the continuum of band states within each bin is 
approximated by a single representative state (the linear com- 
bination of states within the bin that couples to the molecular 
orbitals), Eq. ([T"3l is mapped via a Lanczos transformation to 



Pleads ~ ^ Zj T "(f"°-f r 



»+l.(T 



+ f f ) 



(34) 



representing a semi-infinite, nearest-neighbor tight-binding 
chain to which the impurity couples only at its end site n — 0. 
Since the hopping decays exponentially along the chain as 
t„ ~ DA""/ 2 , the ground state can be obtained via an itera- 
tive procedure in which iteration N involves diagonalization 
of a finite chain spanning sites « < N. At the end of iteration 
N, a pre-determined number of low-lying many-body eigen- 
states is retained to form the basis for iteration N +1, thereby 
allowing reliable access to the low-lying spectrum of chains 
containing tens or even hundreds of sites. See Ref. |53| for 
general details of the NRG procedure. 

For our problem, NRG iteration N - treats a Hamilto- 
nian H = H mo \ + ffmoi-ieads, with N; U2 £ k c eka - in Eq. (O 
replaced by V2/oo-- Since the phonon mode described by H p ^ 
has an infinite-dimensional Hilbert space, we must work in 
a truncated space in which the boson number is restricted to 
n b < N b . 



A. Thermodynamic quantities 

The NRG method can be used to evaluate a thermodynamic 
property X as 



(35) 



where I*!*,,,) is a many-body eigenstate at iteration N having 
energy E m ,p = 1/IcbT, and 



-PE,„ 



(36) 



is the partition function evaluated at the same iteration. For 
a given value of N, Eqs. d35l l and (f36t provide a good 



account^— of X(T) over a range of temperatures around Tn 
defined by k B T N = DA~ N/2 . 

For extensive properties X, it is useful to define the molec- 
ular contribution to the property as 



X 



mol 



x u 



eadsi 



(37) 



where X tot (Xi e ^s) is the total value of X for a system with 
(without) the molecule. In our problem, the local phonon 
mode is treated as part of the host system. Accordingly, we 
define the molecular entropy as 



S mo i(T) - S tot (T) - Si e ads(r) - S p h(T), 



(38) 



where S tol (T) is the total entropy of the system, S\ eai i s (T) is the 
contribution of the leads when isolated from the molecule, and 
S P h{T) is the entropy of the truncated local-phonon system, 
given by 



S P h( T ) = MlnZph(r) - d\nZ ph /8/3], 



with 



,V/. 



(T) = 2 



n b =0 



1 _ e -jS*<uo(iV»+l) 
I — e -jSfiwo 



(39) 



(40) 



Another property of interest is the molecular contribution to 
the static magnetic susceptibility, 

X(T) = J][<*«IS^»>-I<*».|5J*»>I 2 ]^. (41) 

m 

where S z is the total spin z operator, fis is the Bohr magne- 
ton, and g is the Lande g factor (assumed to be the same for 
electrons in the leads and in the molecular orbitals). One can 
interpret |/x mo i| 2 = ?>Tk B x mo \ as the magnitude-squared of the 
molecule's effective magnetic moment. 



B. Linear-response transport properties 

In this paper, we restrict our calculations to equilibrium sit- 
uations in which no external bias is applied. In such cases, 
inelastic transport produced by the e-ph interaction can be 
neglected^ and the linear conductance through the molecule 
can be obtained from a Landauer-type formula 

G(r) = G J {-?£}[-1mT((o,T)]d(o, (42) 



where 



T(oj,T) 



2D 



(43) 
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and Go = 2e 2 /h is the quantum of conductance. The fully 
dressed retarded molecular Green's functions G^iut.T) are 
defined by 

([d^d'liO^e^'dt, (44) 



where (■ ■ ■ ) represents the equilibrium average in the grand 
canonical ensemble and r\ is a positive infinitesimal real num- 
ber. 

As shown for the related problem of two quantum dots con- 
nected in common to a pair of metallic leads^ in the case 
Va = V,c = Vj assumed in the present work, Eq. d42b can be 
recast in the simpler form 

G(T)/G = 7tT c Y j j (~)AZ(a>,T)da>, (45) 

where r c = E„ + Tp and A%(oj, T) = n~ l \m G£(w, T) [defined 
via Eq. ( l44l il is the spectral function for the current-carrying 
linear combination of the a and /? orbitals: 

d m = J] y/fWcdur. (46) 

Within the NRG approach, one can calculate 

A" cc (a>,T) = i 2|(¥ m ,|4j¥ m >| 2 ( e -^» 

X 6 T (io ~{Ert- E m )lh), (47) 

where 6t(co) is a thermally broadened Dirac delta function. 53 
We focus below exclusively on situations where there is no 
magnetic field, and hence A^ c (oj, T) = A cc (oj, T), independent 
of cr. 



IV. NUMERICAL RESULTS 

This section presents NRG results obtained for the Hamil- 
tonian defined by Eqs. (HJ-©, (O, and (O. We take the 
half bandwidth D = 1 as our primary energy scale and adopt 
units in which Ti — kg = g^B — E All results shown 
were obtained for the special case of equal orbital hybridiza- 
tions V a — Vp — V and equal bare Coulomb repulsions 
U a = Up — U' — U . These choices simplify several of the 
effective couplings entering Eq. ([28}, but qualitatively simi- 
lar results are obtained for more general parameter choices. 
The calculations were performed for hybridization V = 0.075 
(hence, E = nV 2 /D 0.0177) and phonon energy ojq = 0.1 
with NRG discretization parameter A = 2.5, allowing up to 
Nb = 60 phonons in the local mode, and retaining 2 000-4 000 
many-body states after each iteration. 

A. Upper orbital far above chemical potential 

We first consider the case of fixed sp = 4 where the up- 
per molecular orbital lies far above the chemical potential 




E 

a 



FIG. 2: (Color online) Variation with orbital energy s a at zero tem- 
perature of (left) the charge (n mol ), and (right) the phonon occupation 
(n b ), for U = (top), U = 5 (middle) and U = 0.5 (bottom). Data 
are for Eg = 4, V = 0.075, A' = A, a> = 0.1, and the four values of 
A 2 /o)o listed in the legend. 

of the leads and therefore contributes little to the low-energy 
physics. This situation, in which the two-orbital model largely 
reduces to the Anderson-Holstein model^^i serves as a ref- 
erence against which to compare cases in which both molec- 
ular orbitals are active. Given that the fi orbital will have neg- 
ligible occupation, the phonon-assisted interorbital tunneling 
i/tun is not expected to play a significant role. In Figs. [2] and 
[3] we show results obtained for A' = A. Figure [4] shows that 
almost identical properties arise for A' = 0. 

Figure [2] shows the total charge of the molecule and the 
phonon occupation as functions of e a , calculated at T = 
for A' = 0, four values of A, and three values of U. First 
consider the case U = of vanishing e-e interactions shown in 
panels (a) and (d). For A — 0, s„ — is a point of degeneracy 
between configurations having molecular charges 0, 1, and 2; 
(n mo i) increases from to 2 over a narrow range As a 4E as 
the a orbital drops below the chemical potential of the leads. 
For A > 0, U = -2A 2 /u>o is negative, and the ground state 
switches from charge to charge 2 around 2s a + U — or 
e a = 2/l 2 /«o- There is a marked decrease in the step width 
As a with increasing A (to be discussed further below). 

It is evident from Figs. |2ja) and |2jd) that changes in the 
ground-state phonon occupation are closely correlated with 
those in the total molecular charge. The prediction of Eq. 
d33b for the case A' - A (hence, A e = 2A and A„ = 0) is 
(rib) = (2/1/wo) 2 («(?)■ Although this relation captures the 
correct trends in the variation of (rib) with s a in Fig. HJd), it 
overestimates the phonon occupation by a significant margin. 
Such deviations are not unexpected, given that Eq. d33l l was 
derived under the assumption that fito^ is the largest energy 
scale in the problem, whereas here sp is the dominant energy 
scale. Empirically, we find that (rib) lies closer to 

rib = (2A/coo) 2 (n e ) 2 , (48) 

which also serves as a lower bound on the phonon occupation. 
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FIG. 3: (Color online) Zero- temperature conductance G vs orbital 
energy e a for (a) U = 0, (b) t/ = 5, and (c) 1/ = 0.5. Data are for 
eg = 4, V = 0.075, /I' = /I, ljq = 0.1, and the four values of A 2 /a)o 
listed in the legend. The inset in (b) shows the data from the main 
panel replotted as G vs e a /U. 



FIG. 4: (Color online) (a) Zero-temperature conductance G vs orbital 
energy e„ for the same parameters as used in Fig. [3jc), except here 
A' = 0. (b) Data for A' = A from Fig.[3lc) replotted as G vs s a /U. (c) 
Data for A' = from (a) replotted as G vs e a /U. 



The error {rib) -fib is largest in the vicinity of the step in (n mo i) 
and vanishes as (« mo i) approaches or 2. For A 2 /o»o = 0.064, 
(n 6 ) - «i < 0.06 both for e a < 0.149 and for s a > 0.151, 
whereas for s a = 0.149945, hi, - 0.88 underestimates (n^,) by 
approximately 0.68. The peak error is smaller for the other 
e-ph couplings shown in Fig. |2td). 

For the case U — 5 of very strong e-e interactions [Fig.[2fb) 
and|2e)], the molecular charge rises from to 1 around the 
point where e a = 0. In contrast with the situation for U — 0, 
the step width As a shows no appreciable change with A. The 
phonon occupation across the charge- 1 is described very well 
by Eq. (g8]>. 

Lastly, Figs. |2]c) and|2jf) show data for U = 0.5, exempli- 
fying moderately strong e-e interactions. With decreasing s a 
(at fixed A), the molecular charge rises in two steps, first from 

to 1 as e a decreases through the chemical potential, and then 
from 1 to 2 as U + s a reaches zero. Just as for U = 5, each 
step has a width As a = 0(T) that is independent of A over 
the range of e-ph couplings shown. The distance along the 
s a axis between the two steps (i.e., the width of the charge- 

1 plateau), which is roughly U — U — 2A 2 /a>o, decreases as 
the e-ph coupling increases in magnitude. Away from the im- 
mediate vicinity of the charge steps, the phonon occupation is 
well-approximated by hb given in Eq. ( l48l i. 

Figure |3]plots the zero-temperature linear conductance as a 
function of s a for the same set of parameters as was used in 
Fig. [2] At T — in zero magnetic field, Eq. ( |45T l reduces to 
G(T = 0)/Go = nr c A cc (Q, 0). In any regime of Fermi-liquid 
behavior, A CC (Q, 0) is expected to obey the Friedel sum-rule, 
implying that n r c A ec (0, 0) = 1 in the wide-band limit where 
all other energy scales in the model are small compared with 
D. This property, which should hold even in the presence of 



e-ph interactions within the molecule, leads to 

G(T = 0) = G sin 2 g<« mol >). (49) 

For (7 = [Fig. [2 a)], we observe a conductance peak at the 
point of degeneracy between molecular charges and 2. This 
is the noninteracting analog of the Coulomb blockade peak 
seen in strongly interacting quantum dots and single-molecule 
junctions above their Kondo temperatures. For A = 0, the 
peak is located at s a = e a — and has a full width As a 
2F, as expected for this exactly solvable single-particle case. 
With increasing A, the conductance peak shifts in position to 
s a - 2A 2 /ol>o (the midpoint of the steps in (n mo i) and («/,)) 
while its width narrows, reflecting the same renormalization 
of r inferred from Figs. [2 a) and[2jd). For all values of A, the 
maximum conductance is G = Go, as predicted by Eq. d49b 
for the point where (n mo t) passes through 1. 

For the interacting cases shown in Figs.[3jb) and[3jc), the 
formation of a many-body Kondo resonance at the chem- 
ical potential leads to a near-unitary conductance at low- 
temperatures T <sz T K over the entire range of s a for which 
("moi) - I- In the case U = 5, no data are shown for 
s a < -0.4, a range in which the Kondo temperature TV is 
so low that the ground-state properties are experimentally in- 
accessible. For both U = 0.5 and U = 5, the conductance 
data show quite good collapse when plotted against s a /U [see 
Fig- 0b) inset and Fig.Ub), respectively]. This scaling con- 
firms that the main features are associated with the conditions 
s a = and e a = -U. 

It has already been remarked that the phonon-assisted in- 
terorbital tunneling is expected to play little role in cases 
where the f3 orbital is far above the Fermi energy. To investi- 
gate this issue, we have compared data for A' — A and A? = 
with all other parameters the same. The molecular charge and 
the phonon occupation are virtually indistinguishable in the 
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two cases. The conductance curves are also very similar, as 
exemplified for U = 0.5 by Figs. [3jc) and IS a). However, 
there are subtle differences that can be highlighted by replot- 
ting G against the scaling variable e a /U. The conductance 
steps in the X = A data [Fig. |4[b)] conform well to the condi- 
tions e a = and e a = —U, but the collapse is even better in 
the case X = [Fig. Etc)]. 

We conclude this discussion of situations in which the /3 
molecular level lies far above the chemical potential by re- 
turning to the narrowing with increasing A at U = of the 
steps in the molecular charge and the phonon occupancy, and 
of the peaks in the linear conductance. This behavior is asso- 
ciated with the presence of a direct crossover of (n mo i) from 
to 2, as opposed to distinct crossovers from to 1, then from 1 
to 2. The same step narrowing can be seen for U > when A 
becomes large enough to suppress the (n mo i) = 1 plateau. (For 
U = 0.5 and X = A, this takes place around A 2 /loo = 0.15, 
considerably larger than any of the values shown in Figs. [2] 
and[3]) 



B. Both orbitals near the chemical potential 

The rich behavior of the model described by Eqs. (HJ-® 
becomes apparent only in the regime where both molecular 
orbitals lie close to the chemical potential and therefore can 
contribute to the low-energy physics. We focus on situations 
in which the orbitals are arranged symmetrically around the 
chemical potential, i.e., Ep = —s a — 6 > 0, a small energy 
scale. (Other arrangements are briefly considered at the end 
of the subsection.) For simplicity, we take X = A throughout. 

We begin by using the transformed Hamiltonian H defined 
in Eq. ( |25T l to calculate the energies of the low-lying states 
of the isolated molecule in the absence of any electron tun- 
neling to/from the leads (i.e., for V = 0). As in the previ- 
ous subsection, we take U a = Up = U' = U, a condition 
that ensures / = W = 0, U± = U\\, and t — -6 in Eq. ( f26l >. 
Then the only explicit e-ph coupling remaining in H appears 
in the term 1 2 ( r(^2^' c <L rC w + H.c). Since \t\ — 6 is small, this 
Hamiltonian term causes only a weak admixture of states with 
nonzero phonon occupancy into the ground state. We there- 
fore focus on the low-lying states of PoH e Po, where H e given 
in Eq. (|26l > is the pure-electronic part of H and Pq projects 
into the zero-phonon sector of the Fock space. Table J] lists 
the low-lying energy eigenstates in this sector for the case 
6 = where both the a and the (3 molecular orbital lie pre- 
cisely at the chemical potential. Also listed are the energies 
of these states including the leading perturbative corrections 
for 5 > 0. These corrections contain a multiplicative factor 
K0|fi ± 2,i<|0>| 2 = exp[-4(/l/w ) 2 ] reflecting the reduction with 
increasing e-ph coupling of the overlap of the phonon ground 
states for Fock-space sectors of different n mo \. 

It can be seen from Table|I]that the singly occupied states of 
lowest energy are l^j 1 ') and \(f>i '), while the doubly occupied 
state of lowest energy is \<j>\ ). One can use the energies of 
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TABLE I: Low-lying eigenstates ]^" mol) > of P H„P where # c is de- 
nned in Eq. J26t > and Pq is a projection operator into the zero-phonon 
sector of the Fock space. Eigenstates for 6 = Sp = —s a = 0, 
U a = Up = U' = U , and A' = A > are grouped according to 
their total electron number n mo |, and specified in terms of operators 
d) y(T defined in Eqs. \2\\ acting on |0>, the vacuum state for elec- 
trons and phonons. Ej™ is the energy of state |0j 1 ) including 
the leading perturbative correction for 6 > 0, expressed in terms of 
U, x = A 2 /it) , andy = a> (6/A) 2 exp[-A(A/a> ) 2 ]. (The nondegener- 
ate perturbation theory used to derive the energies is valid only for 
A * 0.) 

these states to define an effective Coulomb interaction 

U eS = - 2E^ = U-*-* + 5 ^ exp[-4 Wwo ) 2 ]. (50) 
ojq 3A 1 1 1 

In the limit of small 6, this value decreases with increas- 
ing A 2 /(Do at a rate four times greater than the effective 
Coulomb interaction U a [Eq. d!7bH in the a level when the 
j3 level is far above the chemical potential. This property 
of molecules having multiple, nearly-degenerate orbitals en- 
hances the prospects of attaining a regime of negative C/ e ff and 
may have interesting consequences in the area of supercon- 
ductivity. 

Table Q] also indicates that the ground state of the isolated 
molecule crosses from single electron occupancy (for smaller 
\A\) to double occupancy (for larger \A\) at an e-ph coupling de- 
termined by the condition = E^, which yields A 2 /a>u = 
(f//12)[l + 0(5 2 /U 2 )]. As will be discussed in the paragraphs 
that follow, a similar crossing that takes place in the full prob- 
lem with nonvanishing coupling between the molecule and 
the leads results in pronounced changes in the system's low- 
temperature properties. Since the lowest energy of a state hav- 
ing n mo \ > 2 (not shown in Tabled is 3U - 16x-y/12 (which 
exceeds E^ by at least 2U), such states play no role in the 
low-energy physics. 

Figure |3a) shows the evolution of the zero-temperature 
molecular charge (n mo i) with A 2 /u>q for U = 0.5 and four val- 
ues of 6 satisfying 6 <sz U, while Fig. [6] presents the corre- 
sponding occupancies of individual molecular orbitals: (n a ) 
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FIG. 5: (Color online) Variation with A 2 /a>o of (a) the ground-state 
molecular charge (Wm i) = (n e + n a ), (b) the ground-state phonon 
occupation (ni,), (c) the crossover temperature T", and (d) the zero- 
temperature linear conductance G, all calculated for U = 0.5, V = 
0.075, a>o = 0.1, and the four values of 6 = sp = —e a listed in the 




FIG. 6: (Color online) Occupation of individual molecular orbitals 
vs A 2 /u)o for U = 0.5, V = 0.075, too = 0.1, and the four values of 
S = £p = —e a listed in the legend: (a) (n a ) (open symbols) and (np) 
(filled symbols); (b) (n e ) (open symbols) and (n ) (filled symbols). 



and (np) in panel (a), and (n e ) and (n ) in panel (b). We 
begin by considering the special case A. = in which the 
electron and phonon subsystems are entirely decoupled. For 
5 « F = nV 2 /D 0.0177, <« mol > * (n e ) 0.5, which 
may be understood as a consequence of the ground state be- 
ing close to that for U — V = oo and 6 — 0: a product of (1) 
i[cUl r cl/ et - V2c^"J|0> where c m = (2N s T m Zk c ek ^ 
annihilates an electron in the linear combination of left- and 
right-lead states that tunnels into/out of the molecular orbitals, 
and (2) other lead degrees of freedom that are decoupled 
from the molecule. The total charge increases with 5 and ap- 
proaches (rc mo i) = (n a ) = 1 for 6 » T, in which limit the 
large Coulomb repulsion U leads to local-moment formation 



in the a orbital. The local moment is collectively quenched by 
lead electrons, leading to a Kondo singlet ground state. Figure 
|3Jc) shows the characteristic temperature T* of the quenching 
of the molecular spin degree of freedom, determined via the 
standard criterion^ T*x mo i(T*) = 0.0701. This temperature is 
of order F in the mixed-valence limit i «T, but is exponen- 
tially reduced in the local-moment regime 6 » F where T* 
represents the system's Kondo scale. 

Turning on e-ph couplings X = A lowers the energy of the 
even-parity molecular orbital below that of the odd orbital, 
and initially drives the system toward (n e ) = 1, (n ) = 0, 
and toward a many -body singlet ground state formed between 
the leads and a local moment in the even-parity molecular or- 
bital. T* in Fig.|3c) shows an initial decrease with increasing 
A 2 /a>o that is very strong for the smaller values for 6, where 
the e-ph coupling drives the system from mixed valence into 
the Kondo regime. For larger <5, where the system is in the 
Kondo limit even at A = 0, there is a much milder reduction of 
T* caused by the phonon-induced shift of the filled molecular 
orbital further below the chemical potential. 

Upon further increase in the e-ph coupling, {n mo \) and 
T* both show rapid but continuous rises around some value 
A = A x . The crossover value A^/cjq 0.042, which is inde- 
pendent of 6 for 6 <k U, coincides closely with its 6 — > value 
UjYl- 0.0417 for the isolated molecule, where it describes 
the crossing of the singly occupied state \y\ ) and the doubly 
occupied state |i/' ( 1 2) ) (see Table U). For F > 0, the correspond- 
ing energies and £j 2) each acquire a half width F so the 
crossover of the ground-state molecular charge from 1 to 2 is 
smeared over the range \U - 12/t 2 /wo| < 2F. This implies a 
fuU width for the crossover A{A 2 Ioj q ) = 4F/12 = 0.006, an 
estimate in good agreement with Figs. [2 a) and|6] 

In the regime A > A x , minimization of the e-ph energy 
through (n e ) 2, (n„) - (shown in Fig. [6] to hold for all 
the 5 values considered) outweighs the benefits of forming a 
many-body Kondo singlet. Therefore, T* characterizing the 
vanishing of Txmoi ceases to represent the Kondo tempera- 
ture and instead characterizes the scale, of order l2A 2 /a>o - U, 
at which deviations of the molecular occupancies from their 
ground-state values become frozen out. 

Over the entire range of 6 and A 2 /ci>o illustrated in Figs. 
|5]and|6] the ground-state phonon occupation {rib) [Fig-Eb)] 
closely tracks Fit, defined in Eq. d48l ) to within an absolute er- 
ror < (rib) ~ fib - 0.2, an error that peaks around A = A x . 
Similarly, the T — conductance [Fig. [5jd)] is everywhere 
well-described by Eq. ( |49l , reaching the unitary limit Go over 
a window of Kondo behavior for A < A x in which the molecu- 
lar occupancy is 1 , then plunging to zero as the Kondo effect 
is destroyed and the occupancy rises to 2. 

We note that the equilibrium properties shown in Figs. [5] 
and|6]exhibit no special features in the resonant case 5 — 0.05 
in which the molecular level spacing ep - s a exactly matches 
the phonon energy. We expect the resonance condition to play 
a significant role only in driven setups where a nonequilibrium 
phonon distribution serves as a net source or sink of energy for 
the electron subsystem. 

To this point, we have focused entirely on ground-state 
(T = 0) properties. We now briefly discuss the finite- 




FIG. 7: (Color online) Temperature dependence of (a) the molecular 
entropy, (b) temperature times the molecular susceptibility, Tx mo i = 
lA*moil 2 /3, where fi mo \ is the molecule's magnetic moment, and (c) 
the phonon occupation. Data for V = 0.5, S = 0.1, V = 0.075 and 
the four values of A 2 /<±>o listed in the legend. In (a), the horizontal 
dashed lines mark S mo \ = In 2, In 3, and In 5. In (c), the dashed line 
show the occupation of a free phonon mode of energy o>q = 0. 1. 



temperature behavior in situations where the molecular or- 
bitals are arranged symmetrically around the chemical poten- 
tial. Figure |7]plots the temperature dependence of the molec- 
ular entropy, molecular susceptibility, and phonon occupation 
for U = 0.5, 6 = 0.1, V = 0.075, and four different values 
of A 2 /u)q. As long as the temperature exceeds all molecu- 
lar energy scales, the entropy and susceptibility are close to 
the values S mo \ = In 4 and Tx m o\ = 1/8 attained when ev- 
ery one of the 16 molecular configurations has equal occu- 
pation probability, while the phonon occupation is close to 
the Bose-Einstein result for a free boson mode of energy cjq 
[dashed line in Fig. [TJc)]. Once the temperature drops below 
U, most of the molecular configurations (and all with total 
charge n mo \ > 2) become frozen out. For A <K A x (exemplified 
by A 2 /oj = 0.025 in Fig. 13), there is a slight shoulder in the 
entropy around S mo i = In 5 and a minimum in the square of the 
local moment around T^moi = 1/5, the values expected when 
the empty and singly occupied molecular configurations (the 
first five states listed in TableQ]are quasidegenerate. At lower 
temperatures, there is an extended range of local-moment be- 
havior (S mo i = In 2, Txmoi - 1/4) associated with single oc- 
cupancy of the even-parity molecular level (states \<p^ ') and 
\<pi ) )). Finally, the properties cross over below the tempera- 
ture scale T* defined above to those of the Kondo singlet state: 
S mol = 0, TXmol = 0. 

For A just below A x (/l"/wo = 0.0391 in Fig. [7]i there are 
weak shoulders near S mo \ — In 5 and Txmoi = 1/5, as in 
the limit of smaller e-ph couplings. In this case, however, 
these features reflect the near degeneracy of the four « mo i = 1 
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FIG. 8: (Color online) Magnetic moment \r = Txma\ vs scaled tem- 
perature TjT' for U = 0.5, S = 0.1, V = 0.075 and values of A 2 /a> 
spanning the crossover from the Kondo regime to the doubly occu- 
pied regime. The collapse over the range T < 10T" of all curves 
corresponding to A 2 /a> < 0.0391 demonstrates the universal physics 
of the Kondo regime. No such universality is present in the boson- 
dominated limit. 



configurations and the lowest-energy n 



|^j 2) ) in TableU At slightly lower temperatures, the states \<f>P) 
and \<p[ ) become depopulated and the properties drop through 
S mol = In 3 and Txmoi =1/6 before finally falling smoothly 
to zero. Even though there is no extended regime of local- 
moment behavior, the asymptotic approach of S mo \ and Txmoi 
to their ground state values is essentially identical to that for 
A <k A x after rescaling of the temperature by T*. As shown in 
Fig. [8] throughout the regime A < A x , Txmoi follows the same 
function of T/T* for T < 10T*. This is just one manifestation 
of the universality of the Kondo regime, in which Tk = T* 
serves as the sole low-energy scale. 

A small increase in A 2 /cu from 0.0391 to 0.04389, slightly 
above A 2 K /a>o = 0.0417, brings about significant changes in the 
low-temperature properties. While there are still weak fea- 
tures in the entropy at In 5 and ln3, the final approach to the 
ground state is more rapid than for A < A x , as can be seen from 
Fig. [8] Note also the upturn in («/,) as T falls below about 
10J* — a feature absent for A < A x that signals the integral 
role played by phonons in quenching the molecular magnetic 
moment. 

Finally, in the limit A » A x (exemplified by A 2 /(l>q = 0.064 
in Fig. |7J, is by a considerable margin the lowest eigen- 
value of PoH e Pu, so with decreasing temperature 5 mo i and 
Txmo\ quickly approach zero with little sign of any interme- 
diate regime. Even though the quenching of the molecular 
degrees of freedom arises from phonon-induced shifts in the 
molecular levels rather than from a many-body Kondo ef- 
fect involving strong entanglement with the lead degrees, the 
A — > oo ground state is adiabatically connected to that for 
A = 0. 

To conclude this subsection, we consider the effect of shift- 
ing the two molecular levels at a fixed, small energy separation 
£p - s a = 28 through the application of a gate voltage V g that 
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FIG. 9: (Color online) Zero-temperature conductance G as a func- 
tion of s a = -5- V ? (where V g is an applied gate voltage) for the five 
values of A 2 /a> listed in the legend and (a) 6 = 0.025, (b) 6 = 0.05, 
and (c) 6 = 0.1. The other parameters are U = 0.5, V = 0.075, and 
OJn = 0.1. 
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d49l ). A plateau at G - Go spans the range of gate voltages 
within which the total molecular occupancy is (n mo i) ^ 1 
[e.g., compare Figs. [9fc) andfTOlall. while the conductance 
approaches zero for larger V g , where the molecular charge 
vanishes, and for smaller V g , where (« mo i) - 2. 

Once again, we begin by considering the limit A = of 
zero e-ph coupling. For 5 = 0.025 <K F = 0.0177, the steps 
between zero and peak conductance are somewhat broader 
(along the s a axis) than their counterparts in cases where the /3 
molecular orbital lies far above the chemical potential [com- 
pare with Fig. Oc)]. This broadening can be understood as a 
consequence of the step in (n mo i) being split into changes in 
(n a ) and in {np). When 6 » T, the /3 molecular level is essen- 
tially depopulated [as can be seen for the V g = in Fig.[6|a)] 
and the conductance steps narrow to a width similar to that for 
e p = 4. 

Increase of the e-ph coupling from zero results in shifts of 
the occupancy and conductance steps to progressively higher 
values of s a (or to lower values of Vg) that can be attributed 
to the phonon-induce renormalization of the orbital energies 
and of the Coulomb interactions. For 6 = 0.025, the width 
of the (n mo i) =i 1, G - Go plateau is close to the value U e g 
defined in Eq.[50] which approaches U - 8 A 2 /too in the limit 
6 < A 2 /u)q satisfied by the 5 = 0.025 curves in Eq.[9|a). Even 
for the 6 = 0.1 curves shown in Fig. |9jb), the plateau width 
is at least U - 8A 2 /a>o, considerably larger than than its value 
U — U — 2A 2 /cl>o when the orbital /3 lies far above the chem- 
ical potential. The occupancy and conductance plateau might 
be expected to disappear once f/ e ff becomes negative around 
A 2 /loo - U/8 = 0.0625. Indeed, the data for A 2 /lo = 0.064 
in Fig. [9] show a narrow conductance peak that can be asso- 
ciated with the rapid decrease of (n mo i) directly from 2 to 
without any significant range of single occupancy [illustrated 
for <5 = 0.1 in Fig. [Mb)]. 



V. SUMMARY 



FIG. 10: (Color online) (a) Ground-state molecular charge and (b) 
ground-state phonon occupation as functions of s a = -S - V g (where 
V g is an applied gate voltage) for S = 0.1 and the values of A 2 /a>o 
listed in the legend. The other parameters are U = 0.5, V = 0.075, 
and cl)q = 0. 1 . 



causes e, in Eq. ([2]) to be replaced by e, - eV g , and s p in Eq. 
( l26l l to be replaced by s p -eV g . Figure|9]plots the gate-voltage 
dependence of the linear conductance for U = 0.5, five values 
of A' = A, and for 6 = 0.025 [panel (a)], 8 = 0.05 (b), and 
8 — 0.1 (c). Figure [10] shows the corresponding evolution of 
the total molecular charge and the phonon occupation for the 
case 5 = 0.1. In both figures, the quantity plotted along the 
horizontal axis is e a = — S — eV g , which allows direct com- 
parison with the results shown in in Figs.[2jc),[2tf), and[3tc) 
for the regime where the /3 molecular orbital lies far above the 
chemical potential. 

Just as in the other situations considered above, the zero- 
temperature conductance obeys the Fermi-liquid relation Eq. 



We have studied the low-temperature properties of a single- 
molecule junction formed by a two-orbital molecule connect- 
ing metallic leads. We have investigated the low-temperature 
regime of the system using the numerical renormalization 
group to provide a nonperturbative treatment of competing 
strong electron-electron and electron-phonon interactions. In- 
sight into the numerical results has been obtained by consid- 
ering the phonon-renormalization of model parameters iden- 
tified through canonical transformation of the starting Hamil- 
tonian. We have focused on two quite different regions of 
the model's parameter space: (1) In situations where one of 
the two molecular orbitals lies close to the chemical poten- 
tial while the other has a much higher energy, the thermo- 
dynamic properties and linear conductance are very similar 
to those predicted previously for a single-orbital molecule, 
showing phonon-induced shifts in the active molecular level 
and a reduction in the effective Coulomb repulsion between 
electrons on the molecule. (2) In the region where the two or- 
bitals both lie close to the chemical potential, the renormaliza- 
tion of the Coulomb interactions is approximately four times 
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stronger than in the case of one active molecular level, en- 
hancing the likelihood of attaining the interesting regime of 
small or even attractive on-site Coulomb interactions. 

In both regions (1) and (2), electron-phonon interactions 
favor double occupancy of the molecule and are detrimen- 
tal to formation of a molecular local moment and to the 
low-temperature Kondo screening of that moment by elec- 
trons in the leads. With increasing electron-phonon coupling, 
the Kondo effect is progressively destroyed and a phonon- 
dominated nonmagnetic ground state emerges in its place. In 
all the cases presented here, this evolution produces a smooth 
crossover in the ground-state properties. Special situations 
that result in first-order quantum phase transitions between 



Kondo and non-Kondo ground states will be described in a 
subsequent publication. 
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